Stability of fixed points in the (4 + e)-dimensional random field O(A^) spin model for 
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We study the stability of fixed points in the two-loop renormalization group for the random field 
0{N) spin model in 4 + e dimensions. We solve the fixed-point equation in the expansion 
and e expansion. In the large-A*' limit, we study the stability of all fixed points. We solve the 
eigenvalue equation for the infinitesimal deviation from the fixed points under physical conditions 
on the random anisotropy function. We find that the fixed point corresponding to dimensional 
reduction is singly unstable and others are unstable or unphysical. Therefore, one has no choice 
other than dimensional reduction in the large- A'^ limit. The two-loop /3 function enables us to find a 
compact area in the {d, N) plane where the dimensional reduction breaks down. We calculate higher- 
order corrections in the 1/A'' and e expansions to the fixed point. Solving the corrected eigenvalue 
equation nonperturbatively, we find that this fixed point is singly unstable also for sufficiently large 
A'' and the critical exponents show a dimensional reduction. 
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I. INTRODUCTION 

The random field 0{N) spin model is one of the sim- 
plest models with both quenched disorder and spin corre- 
lations. It is a fundamental problem to clarify the critical 
phenomena in this model. Dimensional reduction^ is one 
key to solve this problem. Dimensional reduction claims 
that the critical behavior of the random field 0(iV) spin 
model in d space dimensions is the same as that of the 
pure 0{N) spin model in d — 2 space dimensions. If this 
conjecture is true, all critical exponents of the random 
field spin model in d dimensions should be identical to 
those of the corresponding pure model in d — 2 dimen- 
sions. 

Since several rigorous results for the random field Ising 
model {N — 1 case) indicated the failure of dimensional 
reduction to predict the lower critical dimensionS)2ii4the 
breakdown of dimensional reduction with some approx- 
imation methods was discussed in order to obtain intu- 
itive understanding or quantitative information. Fisher 
calculated the one-loop renormalization group, and he 
pointed out the breakdown of dimensional reduction due 
to the appearance of an infinite number of relevant op- 
erators in 4 -I- e dimensions.^ He showed the existence of 
a fixed point corresponding to the dimensional reduction 
for N > 18, but he argued that this fixed point should 
be unstable as far as the number of spin components, N , 
is finite. Therefore, he concluded that the dimensional 
reduction was not valid near four dimensions. Feldman 



found a nonanalytic fixed point in Fisher's renormaliza- 
tion group for several small He obtained nontrivial 
critical exponents shifted from the predictions of dimen- 
sional reduction. Mezard and Young also suggested the 
breakdown of dimensional reduction by replica symmetry 
breakingji Now, many researchers believe that dimen- 
sional reduction is incorrect in dimensions lower than 6. 

Recently, Tarjus and Tissier studied the critical phe- 
nomena of this model in any dimensions and for any value 
of N by using the nonperturbative renormalization group 
method and the replica methodi^ They show the follow- 
ing relation of the critical exponents of the two-point spin 
correlation function: 



V = V = 



N~2' 



(1) 



predicted by dimensional reduction in a certain region 
in the {d, N) plane. Since this relation seems valid for 
> 18 near four dimensions, the consistency between 
their result and that of the Refs. and|3 should be 
studied. 

To understand the consistency of their works, the de 
Almeida-Thouless criterion^ is applied faithfully to this 
model in a simple 1 /A'^-expansion methodiiS It is shown 
that the saddle point of the auxiliary field is stable in 
the random field 0{N) spin model. Also the stability 
argument by Balents and Fisher for random media^^ is 
applied to Fisher's one-loop renormalization group; then, 
the following two possibilities are indicated liS The fixed 
point corresponding to dimensional reduction is singly 
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unstable, or there is no singly unstable fixed point. Com- 
bining these results leads to relation which is the same 
result obtained by Tarjus and Tissier. 

In this paper, we study the stability of the fixed point 
corresponding to dimensional reduction. Particularly, we 
discuss the physical condition on the deviation from the 
fixed point. Since Fisher did not solve the eigenvalue 
problem completely for the stability around the fixed 
point we solve this problem in a 1/7V expansion. In the 
large- limit, the stability of all fixed points is studied. 
The solution of eigenvalue equations for the deviation 
from the fixed points indicates that the only once unsta- 
ble fixed point shows the critical exponents predicted by 
dimensional reduction. Next, we study this fixed point 
by the two-loop renormalization group obtained by Le 
Doussal and Wiese^- and Tissier and Tarjus. ^'^ The dou- 
ble expansion in 1/iV and e enables us to calculate the 
correction to the fixed point. The solution of the eigen- 
value equation shows that the unstable modes pointed 
out by Fisher are fictitious. Therefore, we conclude that 
the fixed point yielding relation |^ is singly unstable 
for sufficiently large N. This result agrees with that ob- 
tained by Tarjus and Tissies^ and also by the simple 
expansionr^- Furthermore, we calculate higher-order cor- 
rections to the eigenvalues. 

This paper is organized as follows. In Sec. |nl we 
briefiy review Fisher's renormalization group analysis for 
the random field O(A^) spin model in 4 -I- e dimensions.^ 
In Sec. mil we discuss possible singularities of the ran- 
dom anisotropy functions for a physical model. In Sec. 
IIVI we treat the one-loop renormalization group in the 
large- A'^ limit. In this limit, the critical phenomena in 
A + e dimensions are shown to be governed by the fixed 
point which gives the result of dimensional reduction. In 
Sec. we investigate the stability of this fixed point 
corresponding to dimensional reduction by the two-loop 
renormalization group. We show that this fixed point 
is singly unstable on the basis of the physical condition 
on the couphng function discussed in Sec. IIIII Thus, 
we conclude that the prediction of dimensional reduction 
for the critical exponents Q holds for sufficiently large 
N. In Sec. IVll we calculate higher-order corrections to 
the eigenvalue and the exponents of the singularities. In 
Sec. IVIII we summarize our results and discuss some 
problems. 



II. FISHER'S ONE-LOOP RENORMALIZATION 
GROUP ANALYSIS 

In this section, we briefiy review Fisher's argument on 
the instability of the fixed point corresponding to the 
dimensional reduction in the one-loop renormalization 
group for the random field O(A^) spin model in 4 -|- e 
dimensionsi^ 



A. Model 

We consider 0(iV) classical spins S{x) with a fixed- 
length constraint S{x)'^ = 1. To take the average 
over the random field, one introduces replicas S"{x), 
a = 1, . . . ,n. We start from a nonlinear a model of 
the following replica partition function and action: 

/n 
a=l 

2—d r n d 

'' a=l fj.= l 




d'^xY.RiS" ■ S^), (2) 



where a is the ultraviolet cutoff and the parameter T 
denotes the dimensionless temperature. The function 
R{S" ■ S^) represents general anisotropy including the 
random field and all the random anisotropics, and is 
given by 

oo 

where denotes the strength of the random field and 
the /ith rank random anisotropy — 1 is the random 
field, and /i > 2 is the second- and higher-rank ran- 
dom anisotropy). These coupling constants are positive 
semidefinite > 0. 

B. One-loop /3 function 

The (3 function dtR{z) at zero temperature can be ex- 
pressed in the loop expansion 

dtRiz) = l3o[R] + f3i[R] +I32[R] + --- . (4) 

Here, we have defined the scale parameter t which in- 
creases toward the infrared direction. Fisher calculated 
the one-loop /3 function in the following form 

Po[R] = -eR{z), (5) 
l3i[R] = 2{N -2)R'{l)R{z)- {N -l)zR'{l)R'{z) 

+ (1 - z^)R'{l)R"{z) + ^R'{zf{N -2 + z^) 
-R'{z)R"{z)z{l z^) + \R"{zf{l - z^ f{&) 

Expanding R{z) around z = 1, we obtain the one- loop (3 
functions for i?'(l) and 

dtR'{l) = -eR'{l) + {N-2)R'{l)\ (7) 
dtR"{l) = -ei?"(l) + 6i?'(l)i?"(l) 

+{N + l)R"{\f + R'{lf. (8) 
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R"{1) 



' N + 7 




FIG. 1: The renormalization group flow for the couphngs 
R'{1) and R"{1). The fixed point (0, j^yf) is unphysical, be- 
cause the corresponding Ai becomes negative. See Sec. IIV Al 



The /3 functions ijHl and © have two nontrivial fixed 
points (see Fig. QJ: 



iR'il),R'lil)) 

e [N ^ 8) + ^y{N - 2){N - 18) 
^N~2' 2{N - 2){N + 7) ' 

e {N-8)- y/{N -2){N -18) 



(9) 



^N-2' 2{N-2){N + 7) - / ' (^^^ 

The formulas for the critical exponents rj and fj, 

V^R'il), f}^iN-l)R'il)-e, (11) 

enable us to obtain the correlation function critical ex- 
ponents . This result confirms one of the predictions 
by dimensional reduction. From the fixed points © and 
(|10|l . we find that these results are applicable only for 
N > 18. The eigenvalues eAi and eA^ of the scaling 
matrix at the fixed points Q and pO|l are given by 



A± = ± 



TV - 18 



N 



(12) 
(13) 



Thus, the fixed point Q is unstable. The fixed point 
((Tn|l seems to be stable for N > 18. 

However, Fisher claimed an instability of this fixed 
point H1U|I . For the one- loop f3 function in terms 
of differential coefficients it!W(l), (fc = 0,1,2,...), 
he obtained a triangular scaling matrix at the fixed 
point (i?'(l)*,i?':(l)*,...,i?('=-i)(l)*,...), whose diago- 
nal components were given in the following series: 



2fc2 



k{N-l) + 2N- 



N-2 
1-/C+— (fc>3). 



- 1 + kNR'Lil)* 



(14) 



Almost all eigenvalues were indicated to be positive for 
a sufficiently large fc, although one should add a term 



2nkP2Pk missed in Eq. (C6) of his paper.^ Then, Fisher 
concluded that there was no singly unstable fixed point 
and dimensional reduction broke down near four dimen- 
sions. In Sec. however, we show that these infinitely 
many relevant modes are unphysical by solving this eigen- 
value problem completely. In the following sections, we 
carefully reexamine the renormalization group. 



III. ALLOWED SINGULARITIES OF THE 
RANDOM ANISOTROPY FUNCTION 

The fixed-point condition of the renormalization group 
determines properties of the function R{z) . Here we dis- 
cuss possible asymptotic behaviors of R{z) near z = ±1. 
The first derivative of the fixed-point equation with re- 
spect to z is 

[-£ +{N~ 3)R'{l)]R'{z) + zR'{z)^ 
-{N + l)zR'{l)R"{z) + (iV - 3 + Az^)R'iz)R"{z) 
-|-(1 - z^)R'{l)R"'{z) - z{l - z')R'{z)R:"{z) 
-3z(l - z^)R"{z)^ + (1 - zYR"iz)R"'{z) = 0.(15) 

If we assume the asymptotic behavior of i?'(z) near z — I, 

R'iz) ^ R'{1) + C{1- z)'' + ■■■ , (16) 

with < a. To discuss a cuspy behavior of R{z) at 
z = 1, we consider only a < 1. The condition (|15|l gives 
the following constraint: 

[-e+(7V-2)i?'(l)]i?'(l) 

-C2a(4a^ +Aa + N- 1)(1 - z)^"-^ = 0. (17) 
For a 1/2, this constraint gives 
1 



a = -(-l + V2^nV) 



c = o. 



and also 



R'il) = 



N -2 



i?'(l) = 0. 



(18) 
(19) 
(20) 
(21) 



Here, the former case given by Eas. H19|l and 1)20(1 shows 
the dimensional reduction. The formulas for (|ll|l . the 
critical exponents obtained by Feldman,™ enable us to 
obtain the critical exponents 



N -2 



Therefore, no a ^ 1/2 is allowed for any N > 2. For 
a = 1/2, the parameter R'{1) can change continuously 
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depending on the constant C. Therefore, only a = 1/2 
allows divergent R"{1). Only in this case does the non- 
trivial critical behavior differ from the prediction of di- 
mensional reduction. Since the initial value R{z) of the 
renormalization group equation ^ is an analytic func- 
tion, the flow of should diverge for the breakdown 
of dimensional reduction. 

The same discussion for z = —\ can be done. The only 
possible singularity is 

i?'(z) = i?'(-l)-|-C(l-hz)i/2 + ... , 

with C = -R'{1). 

Next, wc consider the renormalization group flow of the 
singular function R'{z). If we assume an initial coupling 



= C(l - z)", 



(22) 



with a > 0, the renormalization group transformation 
generates a term — C^(l — z)'^"^^. If 2q; — 1 < a-namely, 
a < 1-the successive transformations produce less power. 
Eventually, the flow generates a term 

R'{z)^{i~zr' , 

with a' < 0, unless a = i or a > 1. To avoid the flow to 
such unphysical regions, we require a condition 

1 



a = — or a > 1 
2 



(23) 



on the initial function H22|l . Therefore, the allowed func- 
tion has the same singularity as that of the fixed-point 
function. 



IV. LARGE-Af LIMIT 

Here, we take the large- iV limit in the one-loop renor- 
malization group with NR(z) finite and redefine R{z) 
by NR{z) —> R{z). The one-loop /3 function for R{z) 
becomes 



dtR{z) = -eR{z) + 2R'{l)R{z)- zR'{l)R'{z) 
+ \R'{zf. 



(24) 



A. Fixed points 



Following the method given by Balents and Fisher^ii 
we consider the flow equation for R'{z) instead of that 
for R{z). Differentiating the one-loop (3 function with 
respect to z, we have 

dtR'{z) = -eR\z) + R'{l)R'{z)~ zR'{l)R"{z) 

+R'{z)R"{z). (25) 

We redefine the parameters 

R'{z) = eu{z), t' EE et, u{l) = a, (26) 



and we consider the fixed-point equation 

= (a - l)u{z) - zau'{z) + u{z)u'(z). (27) 

Substituting z = 1 into Eq. (|27|l . we have two cases 

a = 0,1, (28) 

for finite u'{l). Solving the differential equation H27|l for 
a = 1, we have two nontrivial solutions: 



R'{z) 



and 



R'{z) = ez. 
The first one indicates 

(Ai,A2) = (6,0), 

(i?'(l),i?"(l)) = (6,0). 



(29) 
(30) 



(31) 
(32) 



Thus, the solution H29|) is the "random field solution." The 
second one indicates 



(Ai,A2) = (0,|), 
(i?'(l),i?"(l)) = (e,6). 



(33) 
(34) 



Thus, the solution l|30|l is not the "random field solu- 
tion" but the "random anisotropy solution." 
For a = 0, the nontrivial solution is 



R'{z) = e(z - 1). 



This indicates 



(Ai,A2) = (-e,- 

(i?'(l),i?"(l)) = (0,6) 



(35) 

(36) 
(37) 



Thus, the solution (|35|l is unphysical. This fact indicates 
that the trivial fixed point i?'(z) = is the unique phys- 
ical solution when a — Q. 
If a 7^ 0,1, 



du{z) (a — 1)m(z) 



dz 



za — u{z) 



(38) 



Taking the inversion, we regard z as a function of u. One 
gets 



dz{u) a z{u) 



1 



du a — I u a— l' 
which is easily integrated. Then, we have 

a/{a-l) 



z = u — (a — 1) 



(39) 



(40) 



Because z{u) takes the maximum value 1 at m = a, m(z) 
is double valued as we show in Fig. |21 It is seen from Eq. 
(|38|l that du/dz is ill defined on u = oz. Therefore the 
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FIG. 2: A schematic graph of u{z). Since the derivative of u 
is ill defined on u = az, the solution terminates on this line. 
The above graph represents two solutions meeting at (1, a). 



lower branch terminates at the origin, so that it should 
be continued to the region — 1 < z < 0. This is possible 
only if a/(a — 1) = n is a positive integer. In this case, 



we have 



> 1. 



Expanding u around a, we have 



u — {a — 1) 



a/(a-l) 



a + (u - a) - (a - 1) 1 + 



a/(a-l) 



1 



2a(a- 1) 



Since — 1 < z < 1, we have 



2a{a~l) 



> 0. 



(41) 



(42) 



Thus, we find that the fixed point a must satisfy a{a — 
1) > 0. Since a < gives a negative exponent e + fj < of 
the disconnected correlation function, a < is excluded 
and we have a > 1. The Schwartz-Soffer inequalitjii^ 
2ri > fj also requires the other constraint a < I. There- 
fore, these nonanalytic solutions (|4(J|I are unphysical. 



B. Stabilities of the fixed-point solutions 



solution is stable, and the positive eigenvalue eA > 
indicates that the fixed-point solution is unstable. Nor- 
malizing v(z) appropriately, we can take v{l) = or 



1. R'{z) = € 

For a — I and u{z) = 1, Eq. 144|) becomes 

b + v'{z){l~z)^ Xv{z), (45) 

where b represents ^(l) taking or 1. When 6 = 0, the 
solution is 



v{z)^C{l-z) 



-A 



(46) 



where A < because of the initial condition b = v{l) = 0. 
On the other hand, when 6 = 1, a general solution is 



( ^ / A-l+c(l-z)-^ 
^(")= ln|l-z 



(A ^ 0), 
(A = 0). 



(47) 



Here the condition 6=1 requires that A = 1 and c = 
0. Since the singularity of the physical deformation v{z) 
satisfies the condition H22|) , the allowed value of A is 



A = --, A<-1 or A = l. 



(48) 



This shows that the fixed-point solution is singly unsta- 
ble. Note that the cuspy deformation from this fixed 
point is irrelevant. 



2. R'{z) = tz 

For a — 1 and u{z) = z, Eq. H44|l becomes 

v{z) = \v{z). (49) 

Then, A = 1, and the fixed point is fully unstable. Note 
that this instability occurs for large N . 



Next, we study the stability of the fixed points. Let 
u{z) be a fixed-point solution: 

Q = u{z){a-l) + u'{z){u{z)-az). (43) 

We consider an infinitesimal deformation u{z) —f u{z) + 
v{z) and a — > a + b, with a finiteness condition 
sup_]^<2<]^ \v{z)\ < GO and b < oo. We study the be- 
havior of the first order in v{z) and b: 

v{z){a — 1) + u{z)b + v'{z)[u{z) — az] 
+u'{z)[viz) - bz] = \v{z). (44) 

eA denotes the eigenvalue of the scaling operator. The 
negative eigenvalue eA < indicates that the fixed-point 



3. R{z) = 

Since a = and u = in this case, Eq. (|44ll is —v{z) = 
Xv{z), which means A = —1 for any v{z); thus the trivial 
fixed point is fully stable. 



V. TWO-LOOP RENORMALIZATION GROUP 
ANALYSIS 

Here, we discuss the stability of the fixed point by solv- 
ing the eigenvalue problem of the scaling operator in the 
two-loop renormalization group. 
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A. The two-loop (5 function 

Recently, Le Doussal and Wiese calculated a two-loop 
(3 function at zero temperature. Independently, Tarjus 
and Tissier obtained the following consistent result 

= \{N-2){{l-z^fR"{zf 

^(1 - z'')[izR'{z) - (2 + z'')R\l)]R" {zf 
-2(1 - z^)[R'{z) ~ zR'{l)]R'{z)R"{z) 
+ (1 - z^)R'{l)R\z)^ + AR'{l)^R{z)} 

- z^m z^)R"'{z) -izR^z) R'{z)Y 

X [-(1 - + zR'{z) - R\l)] 

[(iV + 2)(1 - z^)R"{z) - (37V - 2)zR'{z) 
+8K{N ^2)R{z)], (50) 

where c — limzyi \/l — z^R"{z) and K = 2^a is an un- 
known real number. As discussed in the one-loop f3 func- 
tion in Sec. IIIII the possible singularity in the fixed-point 
function R'{z) ~ (1 — z)" and the possible deformation 
from the fixed point are given by a = ^ or a > 1. Other- 
wise, the function cannot be a fixed point. This condition 
is preserved exactly also in the two-loop f3 function. In 
the following, we present some corrected results based on 
the two-loop renormalization group. 

B. Fixed point in the two-loop /3 function 

First, we discuss a condition on N and d for the exis- 
tence of the fixed point corresponding to the dimensional 
reduction. The two-loop P functions for the two coupling 
constants are 

dtR'il) = -eR'il) + {N -2)R'{lf 

+ (iV-2)i?'(l)3, 
dtR"il) = -eR"{l) + {N + 7)R"ilf + 6R"il)R'{l) 

+R'{1)'^ + 2{5N +l7)R"{lf 

+6{N + 7)R"il)^R'il) 

-6{N - 5)i?"(l)i?'(l)2 -{N- 4)i?'(l)3. 

(51) 

As in the one-loop case, the flow of these coupling con- 
stants does not depend on higher derivative couplings. 
The flow is qualitatively the same as the one-loop case 
depicted in Fig. ^ The fixed point corresponding to the 
dimensional reduction has the following correction 

R'(l) = — ( -^-] . (52) 

The fixed-point equation for i?"(l) becomes a cubic al- 
gebraic equation. For the nonexistence of the solution 



R"{1) in the perturbative region, this cubic equation has 
only one real solution. This condition requires the nega- 
tive discriminant 

for the cubic equation aR"{l)^+bR"{l)'^+cR"{l)+d = 0. 
Expanding the left-hand side of the inequality H53|l in 
TV — 18 and e d — 4, we obtain 

49 

N<18~—e. (54) 
5 

This condition determines the region where the dimen- 
sional reduction breaks. The condition on N for the 
nonexistence of the solution i?"(l) is corrected from 
iV < 18 in the two-loop order. This condition is iden- 
tical to the existence condition of a suitable cuspy fixed 
point obtained by Le Doussal and Wieseji^ which is con- 
sistent also with the phase diagram obtained by Tarjus 
and Tissier. In higher dimensions, dimensional reduction 
occurring is more likely. 

C. Fixed point in the double expansion 

The two-loop (3 function H5U|) enables us to calculate 
the higher-order corrections to the fixed point corre- 
sponding to dimensional reduction. For > 18 — ^e, 
this fixed point may exist, and it can be obtained in a 
double expansion with respect to\/N and e to several or- 
ders. We expand the fixed point up to necessary orders 
to discuss its stability. 

We study the stability of this fixed point in the scaling 
operator. In the two-loop analysis, we define the devia- 
tion function from this fixed point by 

5R'{z) = ^v{z). 

We can calculate the correction of these higher orders to 
the eigenvalue equation for the deviation: 

{l-e){l-zf{l + z)v"{z) 

+ [AT - 4z - 2 + e(2 + 4z)] (1 - z)v' [z) 

+ [2(1 - e)z - NX] v{z) + [N -2 + e{A- z)] v{l) ^ 0. 

(56) 

As discussed in the large- A^ limit, we can determine 
the eigenvalues eA. First, we study the equation for 
w(l) = 0. The solutions of this equation have regular sin- 
gular points z = 1 and —1 for the interval —1 < z < 1. 
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Therefore, we can obtain the solutions in the following 
expansion forms: 

oo 

t;(z) = (l-zr^a„(l-z)" (57) 

n=0 

around z = 1 and 

oo 

z;(z) = (l + 2)'3^6„(l + z)" (58) 

n=0 

around z — —1. Substituting these forms into the eigen- 
value equation l|56() , we require that the coeiEcients of the 
lowest order vanish. This requirement gives the indicial 
equations for the exponents a and /3 

2 (1 - e) - (A^ - 4 + 4e) a + 2 - 2e - iVA = 0, (59) 



which have the solutions 
a = a±{\) 



(60) 



- 4(1 - e) ± ^iV2-8iV(l-e)(l-A) 



4(1-6) 



/3 = - 



N 



2(1 -e) 



or 0. 



(61) 
(62) 



The coefficient of an arbitrary order satisfies the following 
recursion relation: 

fc[2(l - e){k + 2a) ~ N + 4(1 - e)]ofc 
-(1 - e)(a + k){a + k + l)ak-i = 0, 

for fc = 1, 2, 3, . . . . The solution of this recursion relation 
is 

^ (l + a)fc(2 + Q)fc qq 

" (3 + 2a - N/{2 - 2e))k 2^ ' ^ ' 

where the symbol (•••)*: denotes 

{x)k^r{x + k)/T{x). (64) 

The solutions can be expressed in terms of the Gaussian 
hypergeometric function 



v{z, a) 

= {i-zr 

X2-F1 ( 1 + a, 2 + a; 3 + 2a - 



N l~z 



2(1 -e)' 2 



where the generalized hypergeometric function is defined 
by the following series expansion 



iF„ {xi,X2, ■ ■ .,Xm;yi,V2, ■■■,Vn;z) 

{Xl)k{x2)k ■ ■ ■ {X,n)k z'^ 



E 

fc=0 



{yi)k{y2)k - ■ ■{yn)k k\ 



(65) 



The general solution is given by the linear combination 
of the two independent solutions 

v{z) — aQv{z, a+(A)) + Oq w(z, a_(A)), (66) 

where a±{\) is defined by Eq. H61|) . In general, this 
solution has a singularity at z = — 1 corresponding to (3 = 
— 2^Y37y given in Eq. (|^ . To obtain physical solutions, 
we eliminate the singularity at z = — 1 by the suitable 
choice of for a requirement \v{z)\ < 00. Also the 
finiteness of the flow requires a_(A) = ^ or a_(A) > 
1 by the condition (|23|l : then, we obtain the possible 
eigenvalues eA given by 



A 



1 9-9e 
2 

A < -1- 



2N 
8-8e 



N 



(67) 
(68) 



Next we consider the case v{l) 7^ 0. The existence of the 
solution given in the expansion requires vanishing the 
coefficient of f (1), which determines A: 



A = Ai = 1 



e 

N 



(69) 



The corresponding special solution is represented in 
terms of the generalized hypergeometric function 



vi{z) 



2-3e 



2-2e 2-2e 



X3F2 l,l,2,;l-a+(Ai),l-a_(Ai); 



1 - z 



(70) 



The finite solution is given by a linear combination of this 
special solution vi{z) and a solution of the homogeneous 
equation with v{l) = and with the A given in Eq. H69|l : 



viz) = a^viz, a+(Ai)) + Vi{z). 



(71) 



Since the coefficient Oq should be chosen for the cancel- 
lation of the singularity at z = —1, this relevant mode 
is unique. Therefore, the fixed point is singly unstable; 
then, dimensional reduction can be observed in the criti- 
cal exponents 77 and for TV > 18— ^e. This result agrees 
with a nonperturbative renormalization group obtained 
by Tarjus and Tissier^ and a simple 1/N expansionpifi 
All eigenfunctions arc singular at z = 1. It is interesting 
that the eigcnfunction H71I) has an essential singularity 
at 1/N ~ 0. This is because of the fact that we solve the 
eigenvalue equation nonperturbatively without a,l/N ex- 
pansion after the derivation of the eigenvalue equation. 
Note that the two-loop renormalization group is useful to 
show the existence of the relevant mode. The limit e — > 
of the relevant mode H70|l corresponds to the eigcnfunc- 
tion of the one- loop scaling operator with the subleading 
correction in the 1/N expansion. The series expansion 
for the relevant mode seems to be ill defined, since 

N 

lim[l-a+(Ai)] =2-- 

e— >0 2 
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is a negative integer for even N > 6. This apparent 
ill definition disappears by the two-loop or other higher- 
order corrections. 

Here, we comment on the infinitely many relevant 
modes pointed out by Fisher,- To compare our result to 
Fisher's one-loop renormalization group analysis, we take 
the limit e — > in the solution (|61|l . They are included 
in the following series 

Q!+(Afe) = fc— 1, (A; = 3, 4, 5, . . .) and Og" = 0. 

These belong to the eigenvalues eAfc given by 



Afc — 1 — fc 



which are positive for sufficiently large k. These agree 
with the eigenvalues (|14|l obtained by Fisher. Since these 
relevant modes diverge at z = —1, we have eliminated 
them as unphysical modes, as discussed above. 



VI. HIGHER-ORDER CALCULATION 

Here, we calculate higher-order correction in the 1/A'' 
expansion to the results. The fixed point corresponding 
to the dimensional reduction can be obtained in a double 
expansion with respect to 1/N and e to several orders: 

e / 1\ e /I 

R{z) 



N 



+ 



e 



e 

1 



2 

-Z + Z 



1 



z'^ + -z^ - 4z - 



7 . 15 o 35 33 
4 4 4 4 



•(72) 



We study the stability of this fixed point in the scaling 
operator. We can calculate the correction of these higher 
orders to the eigenvalue equation for the deviation: 



--(l-zf(l + z)V"(z) 

l-e+^(3 + 2z) + ^(l + 3z) 

X(l-z)2(l + z)i;"(z) 



(73) 



TV - 4z - 2 + e(2 + 4z) - — (7 + 17z + lOz^) 



+ — (31 + 27z + 18z2)J {:\~z)v\z) 



2(l-e)z+— (6z2 + 8z) 



2N 



(15z2 -I- 21z + 2) - A^A v{z) 



1 



TV - 2 + e(4 - z) - ^(7 + 4z + 3^^) 



e /59 „ iz' 

-TV T + ^^-- 



v(\) = 0. 



As discussed in the previous section, we can determine 
the eigenvalues eA. First we discuss the solution with 
v{\) = 0. To obtain finite solutions in the expansion H57(l 
around z = \ and l|58|) around z = — 1, we should require 
that the coefficients of the lowest order vanish. These 
conditions give the indicial equations for the exponents 
a and ji with two- loop corrections: 



4e 
Tv' 



2 1 -e 



iV 



. . 24 23e\ 

Ar-4 + 4e \ a 

AT iV ' 



14 19e 

_2 - 2e + — - — - iVA = 0, 
TV TV 



4 7e' 

2TV + _ ) /3 = 0. 

TV TV,"^ 



(74) 



The second equation has three solutions for /3 = > 0, 
/? = 0, and f3 — (3- < 0. Since the solution with /?_ 
diverges at z — —1, we should eliminate this solution. 
To this end, we need two independent solutions which 
are finite at z = 1. The existence of two solutions of 



■i or a > 1 requires the eigenvalue eA given by 



1 9 - 9e 
A < -1 



57 - 60e 

2A^2 
48 - 38e 



(75) 
(76) 



Next we consider the case w(l) ^ 0. The existence of the 
solution given in the expansion requires vanishing the 
coefficient of w(l), which determines 



A = 1 



e 

TV 



2e 

TV2 



(77) 



The finite solution is given by a linear combination of 
this special solution of the inhomogeneous equation and 
a solution of the homogeneous equation with v{l) = 
and with the A given in Eq. H77|l . Therefore, the fixed 
point is singly unstable; then, dimensional reduction can 
be observed in the critical exponents r/ and fj for TV > 

18- 



VII. DISCUSSIONS 

In this paper, we have studied the stability of the 
fixed points in the renormalization group for the ran- 
dom field 0(TV) spin model in 4 -I- e dimensions for suf- 
ficiently large TV. We argue physical conditions on the 
random anisotropy function. This argument enables us 
to solve the eigenvalue problem for the infinitesimal de- 
viation from all fixed points. In the large-TV limit, all 
nontrivial fixed points are unstable or unphysical, except 
the fixed point corresponding to dimensional reduction. 
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We also discuss the stability of this fixed point on the 
basis of the two-loop renormalization group obtained by 
Le Doussal and Wiesei^ and Tissier and Tarjusii^ These 
corrected results are essentially the same as those in the 
large- A'^ limit. The fixed point corresponding to dimen- 
sional reduction cannot exist for < 18— ^e. This con- 
dition on (d, N) agrees with the phase diagram obtained 
by Tarjus and Tissier— in a nonperturbative renormal- 
ization group and also with an existence condition of 
the cuspy fixed point obtained by the Le Doussal and 
Wiesei-i^ 

We derive the eigenvalue equation of the scaling op- 
erator in the double expansion in and e on the ba- 
sis of the two-loop renormalization group, and solve it 
nonperturbatively. We show that the unique singly un- 
stable fixed point gives the critical exponents ry and fj 
predicted by dimensional reduction. This result agrees 
with that obtained by Tarjus and Tissier in the nonper- 
turbative renormalization group^ and also with the sta- 
bility of the replica-symmetric saddle-point solution in 
the 1/N expansion.^ - Here, we summarize the physical 
insights obtained from our mathematical results for the 
random field O(A^) spin model in 4 -I- e dimensions at 
zero temperature. There are two phases for sufficiently 
weak randomness. The fully stable fixed point R{z) = 
makes the extended ferromagnetic phase. This phase is 
specified by the nonzero constant order parameter. The 
correlation of the spin fiuctuation obeys the power law 
decay with the same exponent as that in the trivial Gaus- 
sian model. In the disordered phase, the order parameter 
vanishes and the spin correlation becomes short ranged. 
For iV > 18 — ^e, the phase transition between these 



two phases is governed by the singly unstable fixed point 
corresponding to dimensional reduction. Therefore, at 
any phase transition point on the phase boundary, the 
critical exponents predicted by dimensional reduction 
should be observed universally. On the other hand, for 
< 18 — ^e, the fixed point corresponding to dimen- 
sional reduction disappears. In this case, it is believed 
that a cuspy nonanalytic fixed point R'{z) ~ (1 — z)^/^ 
governs the phase transition between the ferromagnetic 
and disordered phases. The breakdown of dimensional 
reduction is observed in the critical exponents— 

It is important to explore a breakdown of dimensional 
reduction in some other observables even for TV > 18 — 
^e. For this problem, the statement on the fixed point 
by Tarjus and Tissier is interesting— They argue that the 
singly unstable fixed point has a weak singularity 

R'{z)^{l~zf/^ (78) 

'At z — 1 for large N . Even if the fixed point corre- 
sponding to dimensional reduction has this singularity, 
we cannot find such an essential singularity in the 1/A'^ 
expansion. Therefore, it is possible that the obtained 
fixed point in the 1/N expansion is not analytic and has 
the weak singularity (|78|l . 
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